Population Genetic Analysis of Paris polyphylla var. yunnanensis Based on cpDNA Fragments.

Paris polyphylla var. yunnanensis is a well-known medicinal plant that is mainly distributed in Southwest China; however, its genetic diversity and biodiversity processes are poorly understood. In this study, the sequences of cpDNA trnL-trnF fragments of 15 wild populations and 17 cultivated populations of P. polyphylla var. yunnanensis were amplified, sequenced, and aligned to study the population genetics of this species. Genetic diversity was analyzed based on nucleotide diversity, haplotype diversity, Watterson diversity, population-level diversity, and species-level genetic diversity. Genetic structure and genetic differentiation were explored using haplotype distribution maps and genetic distance matrices. A total of 15 haplotypes were identified in the 32 populations of P. polyphylla var. yunnanensis. Five unique haplotypes were identified from the fourteen haplotypes of the cultivated populations, while only one unique haplotype was identified from the ten haplotypes of the wild populations. The haplotype richness and genetic diversity of the cultivated populations were higher than those of the wild populations (HT = 0.900 vs. 0.861). In addition, there were no statistically significant correlations between geographic distance and genetic distance in the cultivated populations (r = 0.16, p > 0.05), whereas there was a significant correlation between geographical distance and genetic structure in the wild populations (r = 0.32, p > 0.05), indicating that there was a geographical and genetic connection between the wild populations. There was only 2.5% genetic variation between the wild populations and cultivated populations, indicating no obvious genetic differentiation between the wild and cultivated populations. Overall, the genetic background of the cultivated populations was complex, and it was hypothesized that the unique haplotypes and higher diversity of the cultivated populations were caused by the mixed provenance of the cultivated populations.

Dynamic variation of Paris polyphylla root-associated microbiome assembly with planting years.

MAIN CONCLUSION: P. polyphylla selectively enriches beneficial microorganisms to help their growth. Paris polyphylla (P. polyphylla) is an important perennial plant for Chinese traditional medicine. Uncovering the interaction between P. polyphylla and the related microorganisms would help to utilize and cultivate P. polyphylla. However, studies focusing on P. polyphylla and related microbes are scarce, especially on the assembly mechanisms and dynamics of the P. polyphylla microbiome. High-throughput sequencing of the 16S rRNA genes was implemented to investigate the diversity, community assembly process and molecular ecological network of the bacterial communities in three root compartments (bulk soil, rhizosphere, and root endosphere) across three years. Our results demonstrated that the composition and assembly process of the microbial community in different compartments varied greatly and were strongly affected by planting years. Bacterial diversity was reduced from bulk soils to rhizosphere soils to root endosphere and varied over time. Microorganisms benefit to plants was selectively enriched in P. polyphylla roots as was its core microbiome, including Pseudomonas, Rhizobium, Steroidobacter, Sphingobium and Agrobacterium. The network's complexity and the proportion of stochasticity in the community assembly process increased. Besides, nitrogen metabolism, carbon metabolism, phosphonate and phosphinate metabolism genes in bulk soils increased over time. These findings suggest that P. polyphylla exerts a selective effect to enrich the beneficial microorganisms and proves the sequential increasing selection pressure with P. polyphylla growth. Our work adds to the understanding of the dynamic processes of plant-associated microbial community assembly, guides the selection and application timing of P. polyphylla-associated microbial inoculants and is vital for sustainable agriculture.

Complete genome sequence analysis of a new potyvirus isolated from Paris polyphylla var. yunnanensis.

The complete genome sequence of a new potyvirus from Paris polyphylla var. yunnanensis was determined. Its genomic RNA consists of 9571 nucleotides (nt), excluding the 3'-terminal poly(A) tail, containing the typical open reading frame (ORF) of potyviruses and encoding a putative large polyprotein of 3061 amino acids. The virus shares 54.20%-59.60% nt sequence identity and 51.80%-57.90% amino acid sequence identity with other potyviruses. Proteolytic cleavage sites and conserved motifs of potyviruses were identified in the polyprotein and within individual proteins. Phylogenetic analysis indicated that the virus was most closely related to lily yellow mosaic virus. The results suggest that the virus should be classified as a member of a novel species within the genus Potyvirus, and we have tentatively named this virus "Paris yunnanensis mosaic chlorotic virus" (PyMCV).

Comparative transcriptome analysis reveals key genes for polyphyllin difference in five Paris species.

Paris species accumulate a large amount of steroidal saponins, which have numerous pharmacological activities and have become an essential component in many patented drugs. However, only two among all Paris species. Paris are identified as official sources due to high level of bioactive compounds. To clarify the composition of steroidal saponins and the molecular basis behind the differences between species, we investigated transcriptome and metabolic profiles of leaves and rhizomes in Paris polyphylla var. chinensis (PPC), Paris polyphylla var. yunnanensis (PPY), Paris polyphylla var. stenophylla (PPS), Paris fargesii (PF), and Paris mairei (PM). Phytochemical results displayed that the accumulation of steroidal saponins was tissue- and species-specific. PF and PPS contained more steroidal saponins in leaves than rhizomes, while PPY accumulated more steroidal saponins in rhizomes than leaves. PPC and PM contained similar amounts of steroidal saponins in leaves and rhizomes. Transcriptome analysis illustrated that most differentially expressed genes related to the biosynthesis of steroidal saponins were abundantly expressed in rhizomes than leaves. Meanwhile, more biosynthetic genes had significant correlations with steroidal saponins in rhizomes than in leaves. The result of CCA indicated that ACAT, DXS, DWF1, and CYP90 constrained 97.35% of the variance in bioactive compounds in leaves, whereas CYP72, UGT73, ACAT, and GPPS constrained 98.61% of the variance in phytochemicals in rhizomes. This study provided critical information for enhancing the production of steroidal saponins by biotechnological approaches and methodologies.

RNA-Seq analysis reveals the important co-expressed genes associated with polyphyllin biosynthesis during the developmental stages of Paris polyphylla.

BACKGROUND: Plants synthesize metabolites to adapt to a continuously changing environment. Metabolite biosynthesis often occurs in response to the tissue-specific combinatorial developmental cues that are transcriptionally regulated. Polyphyllins are the major bioactive components in Paris species that demonstrate hemostatic, anti-inflammatory and antitumor effects and have considerable market demands. However, the mechanisms underlying polyphyllin biosynthesis and regulation during plant development have not been fully elucidated. RESULTS: Tissue samples of P. polyphylla var. yunnanensis during the four dominant developmental stages were collected and investigated using high-performance liquid chromatography and RNA sequencing. Polyphyllin concentrations in the different tissues were found to be highly dynamic across developmental stages. Specifically, decreasing trends in polyphyllin concentration were observed in the aerial vegetative tissues, whereas an increasing trend was observed in the rhizomes. Consistent with the aforementioned polyphyllin concentration trends, different patterns of spatiotemporal gene expression in the vegetative tissues were found to be closely related with polyphyllin biosynthesis. Additionally, molecular dissection of the pathway components revealed 137 candidate genes involved in the upstream pathway of polyphyllin backbone biosynthesis. Furthermore, gene co-expression network analysis revealed 74 transcription factor genes and one transporter gene associated with polyphyllin biosynthesis and allocation. CONCLUSIONS: Our findings outline the framework for understanding the biosynthesis and accumulation of polyphyllins during plant development and contribute to future research in elucidating the molecular mechanism underlying polyphyllin regulation and accumulation in P. polyphylla.

Complete genome sequence analysis of Paris alphapartitivirus 1: a novel member of the genus Alphapartitivirus infecting Paris polyphylla var. yunnanensis.

A novel double-stranded RNA (dsRNA) virus, tentatively named "Paris alphapartitivirus 1" (ParAPV1, OL960006-OL960007), was detected in Paris polyphylla var. yunnanensis plants exhibiting leaf chlorosis and shrinkage symptoms in Yunnan. Its complete genome sequence was determined using Illumina and Sanger sequencing. ParAPV1 has a bipartite genome that consists of dsRNA1 (1,917 bp) encoding the viral RNA-dependent RNA polymerase (RdRp), and dsRNA2 (1,818 bp) encoding the putative coat protein (CP). Sequence comparisons showed that the RdRp and CP of ParAPV1 are most similar to those of pear alphapartitivirus (PpPV2), with 69.97% and 54.21% amino acid sequence identities respectively. Phylogenetic analysis of the RdRp amino acid sequences of ParAPV1 and other partitiviruses showed that ParAPV1 cluster with viruses in a clade containing alphapartitiviruses, and that its closest known relatives are PpPV2 (BBA66577) and rose partitivirus (RoPV, ANQ45203S). Taken together, these results suggest that ParAPV1 should be regarded as a new member of genus Alphapartitivirus in the family Partitiviridae. This is the first report of a partitivirus infecting P. polyphylla var. yunnanensis.

Four Novel Mycoviruses from the Hypovirulent Botrytis cinerea SZ-2-3y Isolate from Paris polyphylla: Molecular Characterisation and Mitoviral Sequence Transboundary Entry into Plants.

A hypovirulent SZ-2-3y strain isolated from diseased Paris polyphylla was identified as Botrytis cinerea. Interestingly, SZ-2-3y was coinfected with a mitovirus, two botouliviruses, and a 3074 nt fusarivirus, designated Botrytis cinerea fusarivirus 8 (BcFV8); it shares an 87.2% sequence identity with the previously identified Botrytis cinerea fusarivirus 6 (BcFV6). The full-length 2945 nt genome sequence of the mitovirus, termed Botrytis cinerea mitovirus 10 (BcMV10), shares a 54% sequence identity with Fusarium boothii mitovirus 1 (FbMV1), and clusters with fungus mitoviruses, plant mitoviruses and plant mitochondria; hence BcMV10 is a new Mitoviridae member. The full-length 2759 nt and 2812 nt genome sequences of the other two botouliviruses, named Botrytis cinerea botoulivirus 18 and 19 (BcBoV18 and 19), share a 40% amino acid sequence identity with RNA-dependent RNA polymerase protein (RdRp), and these are new members of the Botoulivirus genus of Botourmiaviridae. Horizontal transmission analysis showed that BcBoV18, BcBoV19 and BcFV8 are not related to hypovirulence, suggesting that BcMV10 may induce hypovirulence. Intriguingly, a partial BcMV10 sequence was detected in cucumber plants inoculated with SZ-2-3y mycelium or pXT1/BcMV10 agrobacterium. In conclusion, we identified a hypovirulent SZ-2-3y fungal strain from P. polyphylla, coinfected with four novel mycoviruses that could serve as potential biocontrol agents. Our findings provide evidence of cross-kingdom mycoviral sequence transmission.

[Phylogeography of Paris poliphylla var. yunnanensis based on chloroplast gene trnL-trnF sequences].

Phylogeography is a research hotspot in the field of the genetic diversity and core germplasm construction of endangered rare plants. Paris polyphylla var. yunnanensis is a rare plant species mainly distributed in China. Wild individuals have been overexploited for the last few decades because of increasing demand for such medicines. Therefore, it is great significance to study the phylogeography of P. poliphylla var. yunnanensis based on chloroplast gene trnL-trnF sequences. In this study, chloroplast genes trnL-trnF were used in the phylogeography analysis of 15 wild and 17 cultivated populations of P. polyphylla var. yunnanensis. This study revealed that based on the results of neutrality tests and mismatch analysis, the rapid expansion of wild population has not been detected in P. polyphylla var. yunnanensis. After aligning and sorting the obtained cpDNA sequences, a total of 15 haplotypes were detected in all 32 populations. One haplotype was unique to the wild population, and 5 haplotypes were unique to the cultivated population. It can be seen that the haplotype richness of cultivated population was higher than that of wild population. The wild populations of P. polyphylla var. yunnanensis were divided into two groups according to evolutionary relationship of haplotypes and distribution map of haplotypes. The haplotype of branch Ⅰ was mainly distributed in Guizhou, and the haplotype of branch Ⅱ was located in Yunnan and Huidong, Sichuan. Therefore, it's speculated that Guizhou and the west Yunnan region may be glacial refuge in the evolutionary history of wild populations of P. polyphylla var. yunnanensis, and in order to protect the wild resources more effectively, wild populations of P. polyphylla var. yunnanensis in these two areas should be included in the protection zone.

A cycloartenol synthase from the steroidal saponin biosynthesis pathway of Paris polyphylla.

Steroidal saponins named polyphyllin are the major active components of Paris polyphylla. Cycloartenol synthase (CAS) is a key enzyme that catalyzes the formation of the sterol scaffold. In this study, we cloned a putative CAS gene from Paris polyphylla. Heterologous expression in yeast indicated that PpCAS can convert 2,3-oxidosqualene into cycloartenol. qRT-PCR analysis showed that the expression of PpCAS was highest in leaves and lowest in roots. To our best knowledge, this is the first report of the functional characterization of cycloartenol synthase from Paris polyphylla, which lays the foundation for further analysis of the biosynthesis pathway of polyphyllins.[Formula: see text].

Transcriptome analyses of Paris polyphylla var. chinensis, Ypsilandra thibetica, and Polygonatum kingianum characterize their steroidal saponin biosynthesis pathway.

Steroidal saponins, one of the most diverse groups of plant-derived natural products, elicit biological and pharmacological activities; however, the genes involved in their biosynthesis and the corresponding biosynthetic pathway in monocotyledon plants remain unclear. This study aimed to identify genes involved in the biosynthesis of steroidal saponins by performing a comparative analysis among transcriptomes of Paris polyphylla var. chinensis (PPC), Ypsilandra thibetica (YT), and Polygonatum kingianum (PK). De novo transcriptome assemblies generated 57,537, 140,420, and 151,773 unigenes from PPC, YT, and PK, respectively, of which 56.54, 47.81, and 44.30% were successfully annotated, respectively. Among the transcriptomes for PPC, YT, and PK, we identified 194, 169, and 131; 17, 14, and 26; and, 80, 122, and 113 unigenes corresponding to terpenoid backbone biosynthesis; sesquiterpenoid and triterpenoid biosynthesis; and, steroid biosynthesis pathways, respectively. These genes are putatively involved in the biosynthesis of cholesterol that is the primary precursor of steroidal saponins. Phylogenetic analyses indicated that lanosterol synthase may be exclusive to dicotyledon plant species, and the cytochrome P450 unigenes were closely related to clusters CYP90B1 and CYP734A1, which are UDP-glycosyltransferases unigenes homologous with the UGT73 family. Thus, unigenes of ß-glucosidase may be candidate genes for catalysis of later period modifications of the steroidal saponin skeleton. Our data provide evidence to support the hypothesis that monocotyledons biosynthesize steroidal saponins from cholesterol via the cycloartenol pathway.

Transcriptome-Wide Identification and Prediction of miRNAs and Their Targets in Paris polyphylla var. yunnanensis by High-Throughput Sequencing Analysis.

Long dormancy period of seeds limits the large-scale artificial cultivation of the scarce Paris polyphylla var. yunnanensis, an important traditional Chinese medicine. Characterizing miRNAs and their targets is crucial to understanding the role of miRNAs during seed dormancy in this species. Considering the limited genome information of this species, we first sequenced and assembled the transcriptome data of dormant seeds and their seed coats as the reference genome. A total of 146,671 unigenes with an average length of 923 bp were identified and showed functional diversity based on different annotation methods. Two small RNA libraries from respective seeds and seed coats were sequenced and the combining data indicates that 263 conserved miRNAs belonging to at least 83 families and 768 novel miRNAs in 1174 transcripts were found. The annotations of the predicted putative targets of miRNAs suggest that these miRNAs were mainly involved in the cell, metabolism and genetic information processing by direct and indirect regulation patterns in dormant seeds of P. polyphylla var. yunnanensis. Therefore, we provide the first known miRNA profiles and their targets, which will assist with further study of the molecular mechanism of seed dormancy in P. polyphylla var. yunnanensis.

Hybridization and reproductive isolation between diploid Erythronium mesochoreum and its tetraploid congener E. albidum (Liliaceae).

Polyploidy has played an important role in angiosperm diversification, but how polyploidy contributes to reproductive isolation remains poorly understood. Most work has focused on postzygotic reproductive barriers, and the influence of ploidy differences on prezygotic barriers is understudied. To address these gaps, we quantified hybrid occurrence, interspecific self-compatibility differences, and the contributions of multiple pre- and postzygotic barriers to reproductive isolation between diploid Erythronium mesochoreum (Liliaceae) and its tetraploid congener Erythronium albidum. Reproductive isolation between the study species was nearly complete, and naturally occurring hybrids were infrequent and largely sterile. Although postzygotic barriers effected substantial reproductive isolation when considered in isolation, the study species' spatial distributions and pollinator assemblages overlapped little, such that interspecific pollen transfer is likely uncommon. We did not find evidence that E. albidum and E. mesochoreum differed in mating systems, indicating that self-incompatibility release may not have fostered speciation in this system. Ultimately, we demonstrate that E. albidum and E. mesochoreum are reproductively isolated by multiple, hierarchically-operating barriers, and we add to the currently limited number of studies demonstrating that early acting barriers such as pollinator-mediated isolation can be important for effecting and sustaining reproductive isolation in diploid-polyploid systems.

Conditions in home and transplant soils have differential effects on the performance of diploid and allotetraploid anthericum species.

Due to increased levels of heterozygosity, polyploids are expected to have a greater ability to adapt to different environments than their diploid ancestors. While this theoretical pattern has been suggested repeatedly, studies comparing adaptability to changing conditions in diploids and polyploids are rare. The aim of the study was to determine the importance of environmental conditions of origin as well as target conditions on performance of two Anthericum species, allotetraploid A. liliago and diploid A. ramosum and to explore whether the two species differ in the ability to adapt to these environmental conditions. Specifically, we performed a common garden experiment using soil from 6 localities within the species' natural range, and we simulated the forest and open environments in which they might occur. We compared the performance of diploid A. ramosum and allotetraploid A. liliago originating from different locations in the different soils. The performance of the two species was not affected by simulated shading but differed strongly between the different target soils. Growth of the tetraploids was not affected by the origin of the plants. In contrast, diploids from the most nutrient poor soil performed best in the richest soil, indicating that diploids from deprived environments have an increased ability to acquire nutrients when available. They are thus able to profit from transfer to novel nutrient rich environments. Therefore, the results of the study did not support the general expectation that the polyploids should have a greater ability than the diploids to adapt to a wide range of conditions. In contrast, the results are in line with the observation that diploids occupy a wider range of environments than the allotetraploids in our system.

[Phenotypic Trait Variation, Principal Component, Correlation and Path Analysis of Paris polyphylla var. yunnanensis].

OBJECTIVE: To investigate the phenotypic trait variation range of Paris polyphylla var. yunnanensis, and to look for phenotypic traits closely related with yield, in order to provide reference for its breeding in the future. METHODS: Wild plants of Paris polyphylla var. yunnanensis populations cultivated for three years in Xishuangbanna, Yunnan Province were chosen. The plants of Paris polyphylla var. yunnanensis were randomly labelled, its phenotypic traits were observed and analyzed by variation, principal component, correlation, variance and path analysis. RESULTS: Eleven phenotypic traits variation of the populations of Paris polyphylla var. yunnanensis were great, mainly distributed in yield, growth and genetic characteristics. Phenotypic trait of different growth years (4 - 10 years) including stem height, stem diameter, flower stem height, flower stem diameter, leaf number, leaf length and leaf width had no significant difference(P > 0. 05). The underground phenotypic traits including root weight, root length and root diameter increased significant (P <0. 05) with plant growth years. Effect of the other phenotypic traits on root weight was as follows: root diameter (0. 2550) > root growth years (0. 1183 ) > root length(0. 0202) > stem diameter(0. 0081) > stem height (0. 0044) > leaf number (0. 0025) > leaf length (0. 0018) > leaf width (0. 0004) > flower stem height (0. 0003) > flower stem diameter (0. 0001). CONCLUSION: The phenotypic traits on yield, growth and genetic characteristics have rich variations, and the phenotypic traits on the ground have no relationship with plant growth years, the plant growth years mainly affect the underground phenotypic traits. When choosing different root types, root diameter can be used as the preferred character, stem diameter as the second preferred traits when the root traits information not be easily obtained.

Microbial diversity in Paris polyphylla var. yunnanensis rhizomes of varying ages.

Endophyte microorganisms live inside plants without causing them any apparent damage. Recently, endophytic microorganisms have attracted attention because they can produce bioactive compounds of biotechnological interest. The endophytic microorganisms in Paris polyphylla var. yunnanensis (Liliaceae) - a species used since antiquity in traditional Chinese medicine - are under scrutiny because they may be responsible for producing the bioactive metabolites associated with the plant. The levels of bioactive metabolites in the rhizomes of P. polyphylla increase with rhizome age. To elucidate the roles played by endophytes in the accumulation of bioactive metabolites, we investigated the community structure and diversity of the endophytic microorganisms in P. polyphylla rhizomes of different ages (4, 6, and 8 years) using 16S rRNA and internal transcribed spacer (ITS) sequence analysis. 16S rDNA amplicon pyrosequencing revealed that the number of operational taxonomic units was lower in the 8-year-old samples than in the other samples. A total of 28 phyla were observed in the P. polyphylla samples and the predominant bacteria were of the Cyanobacteria and Proteobacteria phyla. Moreover, the percentage of Cyanobacteria increased with rhizome age. Similarly, ITS1 amplicon pyrosequencing identified developmental changes in the most abundant fungal classes; some classes were more prevalent in the 8-year-old rhizomes than in younger rhizomes, indicating the importance in secondary metabolism in older rhizomes. Our study showed that endophyte microorganism diversity and prevalence depend on P. polyphylla rhizome age. There was also an indication that some endophyte microorganisms contribute to the higher saponin content in older P. polyphylla specimens.

[Construction and preliminary analysis of a full-length cDNA library for Paris polyphylla var. yunnanensis].

OBJECTIVE: A full-length cDNA library of Paris polyphylla var. yunnanensis was constructed in order to research the genes relating to growing development and the genes regulation of its secondary metabolite biosynthesis. METHODS: The total RNA was extracted from Paris polyphylla var. yunnanensis using modified Trizol method. The SMART (switching mechanism at 5' end of RNA transcript )technology was appliedl to construct the full-length cDNA library. The library titer,recombinant rate and length of insert fragments were determined,the sequences of the library were analyzed by Blastx and were compared to GenBank database. RESULTS: The capacity of the library was 2. 5 x 107 cfu/mL, the recombinant rate was 98.5% and the average size of the inserted fragment was 1.5 kb. 9 ESTs (Expressed Sequence Tags) were relating to growing development and 5 ESTs were relating to regulation of secondary metabolite biosynthesis among 149 ESTs obtained from 192 clones sequenced. CONCLUSION: A full-length cDNA library of Paris polyphylla var. yunnanensis is constructed by SMART technology successfully, and the library has enough capacity, high recombinant rate and long insert fragment for the further research to screen and identify the functional genes of Paris polyphylla var. yunnanensis.

Evolutionary relationships in the medicinally important genus Fritillaria L. (Liliaceae).

Fritillaria (Liliaceae) is a genus of approximately 140 species of bulbous perennial plants that includes taxa of both horticultural and medicinal importance. As well as being commercially valuable, Fritillaria species have attracted attention because of their exceptionally large genome sizes, with all values recorded to date in excess of 30Gb. Despite such interest in the genus, phylogenetic relationships between the majority of species have remained untested. Here we present the first phylogenetic reconstruction of relationships to encompass most of the currently recognised species diversity in the genus. Three regions of the plastid genome were sequenced in 117 individuals of Fritillaria, representing 92 species (c. 66% of the genus) and in representatives of nine other genera of Liliaceae. Eleven low-copy nuclear gene regions were also screened in selected species for their potential utility. Phylogenetic analysis of a combined plastid dataset using maximum parsimony and Bayesian inference provided support for the monophyly of the majority of currently recognised subgenera. However, subgenus Fritillaria, which is by far the largest of the subgenera and includes the most important species used in traditional Chinese medicine, is found to be polyphyletic. Moreover, several taxa that were represented by multiple individuals show evidence of species non-monophyly. The Japanese endemic subgenus Japonica, which contains the species with the largest recorded genome size for any diploid plant, is resolved as sister to the predominantly Middle Eastern and Central Asian subgenus Rhinopetalum. Whilst relationships between most of the major Fritillaria lineages can now be resolved, our results also highlight the need for data from additional independently evolving loci; an endeavour that may be particularly challenging in light of the huge nuclear genomes found in these plants.

Origin and diversification of the Milla Clade (Brodiaeoideae, Asparagaceae): a Neotropical group of six geophytic genera.

The Milla clade currently comprises six genera of geophytic plants distributed from Arizona to Guatemala. Three genera (Behria, Jaimehintonia and Petronymphe) are monotypic while the remaining genera (Bessera, Dandya and Milla) contain from two to ten (Milla) species. Parsimony, Maximum Likelihood and Bayesian Inference analyses were conducted with plastid and nuclear DNA sequences from a total of 181 plants belonging to 15 species in all six genera. Molecular dating was performed under a relaxed clock model. We examined the phylogenetic relationships of the genera and species, estimated origin-divergence times for the clade and genera and determined the ancestral distribution area of the clade by optimizing ancestral areas given current biogeographic distributions. The phylogenetic results suggest that final decisions on limits of the six genera in the Milla clade will have to be established until further taxonomic work is completed for Milla, in particular for the group of populations included under the name M. biflora. The later genus is rendered polyphyletic by other genera of the family. The origin of the Milla clade is estimated at 15.8Ma. Ancestral area of the clade most likely was located in the California Floristic Province and dispersal occurred most likely to the Chihuahuan-Coahuila Plateaus and the Trans-Mexican Volcanic Belt and from there to Baja California and the Sierra Madre del Sur. Two hypotheses that need further testing are proposed to explain complex relationships of genera and polyphyly of Milla, one in relation to fragmentation of populations and pollinator shifts and another suggesting that populations remained in refugia in the Trans-Mexican Volcanic Belt.

Phylogenetic reappraisal of Allium subgenus Cyathophora (Amaryllidaceae) and related taxa, with a proposal of two new sections.

The phylogeny of subgenus Cyathophora and representatives of its closely related taxa within Allium were reconstructed based on nrDNA ITS and two plastid fragments (trnL-F and rpl32-trnL). The constructed phylogenies indicated that subgenus Cyathophora was not monophyletic and to be split in three parts positioned in different clusters. Allium kingdonii was unequivocally placed within subgenus Amerallium and formed an immediate sister relationship with New World Amerallium clade, suggesting an unexpected intercontinental disjunct distribution. For another, Allium trifurcatum was firmly nested within subgenus Butomissa next to A. tuberosum and A. ramosum, but it is distinctly different morphologically from the latter by thinly leathery bulb tunics, uniovulate locule and obviously 3-cleft stigma. Based on the geographic features, morphological and molecular evidences, two new sections, Kingdonia X.J.He et D.Q.Huang for A. kingdonii and Trifurcatum X.J.He et D.Q.Huang for A. trifurcatum, were proposed. The remaining three species of subgenus Cyathophora formed a well-defined clade, and the phylogenetic relationships among them recovered were consistent with previous findings. In addition, A. weschniakowii and A. subtilissimum were proven to be a member of subgenera Rhizirideum sensu stricto (s. str.) and Cepa, respectively, rather than subgenera Cepa and Polyprason previously proposed. Section Rhizomatosa represented by A. caespitosum should be subsumed within section Caespitosoprason of subgenus Rhizirideum s. str.

De Novo transcriptome sequencing reveals important molecular networks and metabolic pathways of the plant, Chlorophytum borivilianum.

Chlorophytum borivilianum, an endangered medicinal plant species is highly recognized for its aphrodisiac properties provided by saponins present in the plant. The transcriptome information of this species is limited and only few hundred expressed sequence tags (ESTs) are available in the public databases. To gain molecular insight of this plant, high throughput transcriptome sequencing of leaf RNA was carried out using Illumina's HiSeq 2000 sequencing platform. A total of 22,161,444 single end reads were retrieved after quality filtering. Available (e.g., De-Bruijn/Eulerian graph) and in-house developed bioinformatics tools were used for assembly and annotation of transcriptome. A total of 101,141 assembled transcripts were obtained, with coverage size of 22.42 Mb and average length of 221 bp. Guanine-cytosine (GC) content was found to be 44%. Bioinformatics analysis, using non-redundant proteins, gene ontology (GO), enzyme commission (EC) and kyoto encyclopedia of genes and genomes (KEGG) databases, extracted all the known enzymes involved in saponin and flavonoid biosynthesis. Few genes of the alkaloid biosynthesis, along with anticancer and plant defense genes, were also discovered. Additionally, several cytochrome P450 (CYP450) and glycosyltransferase unique sequences were also found. We identified simple sequence repeat motifs in transcripts with an abundance of di-nucleotide simple sequence repeat (SSR; 43.1%) markers. Large scale expression profiling through Reads per Kilobase per Million mapped reads (RPKM) showed major genes involved in different metabolic pathways of the plant. Genes, expressed sequence tags (ESTs) and unique sequences from this study provide an important resource for the scientific community, interested in the molecular genetics and functional genomics of C. borivilianum.